Second virial coefficient for a d-dimensional Lennard-Jones (2n-n) system 
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This note examines the second virial coefficient for an imperfect 
gas subject to a 2n-n interparticle potential in any dimension d 
between and n. A compact analytic expression is presented for 
this quantity which shows that, apart from a numerical factor, its 
temperature dependence is a universal function parameterized by 
d/n. 
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The virial expansion has played an important role in studying imperfect 
gases and has recently been invoked for the examination of systems, such as 
alkanes absorbed at interfaces [1], in two dimensions. It is therefore of some 
interest to look at the dimensional dependence of the coefficients. In this 
note an analytic expression is presented for the second virial coefficient in 
the case of the Lennard- Jones 2n-n potential 

v(r) = nC-? n - ftl (1) 

for spatial dimension < d < n, where it is defined. 
The second virial coefficient is given by 

B d (T) = X - J d d r[l - exp(-V(r)/kT)\. (2) 

Since V(r) is spherically symmentric, the angular integrals in hyperspherical 
coordinates can be carried out to give 

B ^=Tmr rd ~ lll - e ~ v(T)/lT]dr - (3) 

After scaling out a, integration by parts, and an elementary change of inte- 
gration variable, (3) reduces to 

B d (T) = C d [2h_ d/6 - J_ d/6 ] (4) 

where 

a = Ae/kT 

Cd = mdjY) (5) 



I u = a [°° x v e~ a{ - x2 - x) dx. 
Jo 



The latter integral does not appear to have been tabulated, but is easly eval- 
uated as follows: Replace exp(ax) by its Taylor series and integrate term-by- 
term using Euler's gamma function integral to get an infinite series summed 
with respect to n, say. Then break this up into even (n = 2k) and odd 
(n = 2k + 1) terms. In these series replace n\ by A k k\(l/2) k and 4 k k\(3/2) k , 
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respectively. The resulting series are seen to be Kummer functions [2] and we 
obtain 

v^ + l)M(^ + l,|^)]. (6) 

Therefore, for < d < n, in terms of the dimensionless temperature r = 
A/a = kT/e, and p = d/2n 

d„d/2od/ra 1 

^(r) = -^ w r- 1 r(-p)M(-p,-,0 + 

2r- 1/2 r(^- P)M{\- p\t-*)]= (7) 



a d ir 



d / 2 2 2p 1 1 



T ->e 1 ' r \r(-p)M{p + ^,-l/T)+ 



nT(d/2) vr 2' 2 

A r( l_ p)M (p+l,|-l/r)] 
(7 <i 7r (d-i)/2 2/ , i 



In deriving (7) the identities 

. ,/ d 1 . ^ , d 3 s , d 1 . 

M l , -, z) - 2zM 1 , -, z) = M , -, z) 

v 2n' 2' ; v 2n' 2' ! v 2n' 2' 7 

w 3 d 3 . ,,.1 d 1 . d 3 . 

1 - - M - - — , -, z) - M - - — , -, z) = M - - — - z) 8 

v n ; v 2 2n 2 ; v 2 2n 2 n v 2 2n 2 ; w 

and Kummer's first identity have been used. The last expression in (7) is 
particularly attractive, but as indicated below, the parabolic cylinder func- 
tion can be problematic and the expressions using Kummer functions are 
probably better for numerical work. From Eq.(7) we have the interesting 
fact that, apart from a constant prefactor, the scaled temperature depen- 
dence of the second virial coefficient is a universal function of d/n. Thus, the 
Boyle temperature ( at which the second virial coefficient vanishes), is given 
by the zeros of the parabolic cylinder function. Listed below are the Boyle 
temperatures for the L-J 12-6 potential in various dimensions: 
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d r B d r B 
0.5 0.45812 3.5 5.35369 

1 0.72701 4 9.01642 
1.5 1.07421 4.5 17.1531 

2 1.56031 5 41.0597 
2.5 2.28098 5.5 173.8854 

3 3.41793 

The Boyle temperature is therefore also universal with respect to d/n, so this 
table lists the Boyle temperatures for d/n = 1/12, . . . 11/12. 

Finally we present the Joule-Thompson coefficient Jd,n(j) = rdBd-nij) / 'dr- 
B d , n (r). 

aV ,2 2d ,2n d d d i 

Jd ' n{T) = 6nT(d/2) T {3r( -2^ )[(1 + 2^ )M( ^'2' r } " 

~M(l-^ It -!)] + 
t 2n In 1 

2t- 1 / 2 F(--—)[3(-+—)M{--— ),^, T - 1 )+-(--— )M(-- — A,^ 1 )]}. 
v 2 2n n v 2 2n V 2 2n h 2' ; t2 2n V 2 2n' 2' ;JJ 

( _ (12) 

A similar calculation for the case d = 3 was carried out a number of years 
ago by Garrett [4] who expressed his result in terms of a sum of parabolic 
cylinder functions D v (—z). This function appears to be an orphan among 
special functions, in that there is very little information about it in standard 
references such as [3]. What there is concerns the function D u (z). In spite of 
the similarity in notation, the two are independent solutions of the parabolic 
cylinder equation and their relationship is not a simple one. For example, 
for real z, 

D 1/2 (z) = -L(z/2f 2 [K l4 (z 2 /4) +Ks(z 2 /A)} (13) 

while 

D 1/2 (-z) = ^"[Uy/A) + 7_ f (z 2 /4) - hy/A) - I_y/A)\. (14) 

On the other hand, the Kummer function used in (7) (also denoted M(a, b, c) = 
iFi(a; b; c)) is standard in mathematical packages and does not suffer from 
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this problem. For d — 3 (7) is consistent with Garrett's formula. Recently, 
Vargas et al.[5] have treated B 3t6 by an improvement of Garrett's approach 
and have expressed the virial coefficient in terms of modified Bessel functions 
of fractional order. By making use of (13) and (14), (7) can be transformed 
into their expression. Finally, the method used here for the 2n — n potential: 
series expansion, decomposition of the sum into relatively prime summands 
and the use of the Gauss multiplication formula, can be applied to the general 
m—n potential and leads to a closed form expression as a sum of lcm(m, n)/n 
hypergeometric functions. Thus, for the 12 — 5 potential considered in [1], 
the second virial coefficient was expressed as a sum of 12 Kummer functions, 
which proved to be of little practical value. 
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